Dynamics of Surfactant Driven Fracture of Particle Rafts 
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We investigate the dynamic fracture of a close-packed monolayer of particles, or particle raft, 
floating at a liquid-gas interface induced by the localised addition of surfactant. Unusually for 
a two-dimensional solid, our experiments show that the speed of crack propagation here is not 
affected by the elastic properties of the raft. Instead it is controlled by the rate at which surfactant 
is advected to the crack tip by means of the induced Marangoni flows. Further, the velocity of 
propagation is not constant in time and the length of the crack scales as £ 3//4 . More broadly, this 
surfactant induced rupture of interfacial rafts suggests ways to manipulate them for applications. 



The curious behavior of particulate interfaces that sep- 
arate liquids is a source of interesting questions at the in- 
tersection of hydrodynamics and elasticity. For instance, 
liquid drops coated with a fine hydrophobic powder be- 
come non- wetting 1], forming an artificial analog of a 
much older solution stumbled upon by insects 2]. Sim- 
ilarly, the addition of particles to the surface of liquid 
drops prior to coalescence stabilises the coalesced drops 
to the common pinch-off instability and can lead to re- 
versible morphological instabilities such as buckling when 
subject to pressure 0. Particle covered interfaces also 
occur at an intermediate stage during the production of 
colloidosomes 4], armored bubbles p and porous par- 
ticle aerosols for drug delivery From both a funda- 
mental and technological point of view these particulate 
interfaces pose a number of problems and opportunities. 
Understanding their rheological properties has the po- 
tential to open up ways of controlling them, but remains 
mostly an open question. 

Recent experiments 0, Q show that a densely packed 
monolayer of particles at an interface (a 'particle raft') 
has many of the characteristics of a two-dimensional lin- 
ear elastic solid in certain regimes. For example, un- 
der compressive loading, a particle raft statically buck- 
les out of the plane demonstrating that collectively its 
constituent particles possess a non-zero shear modulus 
G ~ <y/d 8], with 7 the surface tension coefficient of 
the pure liquid-gas interface and d the particle diameter. 
This ability to sustain finite shear stresses arises from a 
combination of capillary forces and the short range steric 
constraints due to particle-particle contact, and is seen 
in a variety of similar systems such as armored bubbles, 
drying colloidal drops and suspensions 0, EJ [HI E| > al- 
though understanding the kinetics of onset of this elastic 
behavior constitutes work in progress. Going beyond the 
linear elastic behavior, the ability to sustain finite shear 
stresses suggests that these rafts should also be able to 
sustain fractures which relieve stresses primarily in one 




FIG. 1: Time series showing the branching fracture of an in- 
terfacial particle raft consisting of 100/xra Pliolite particles 
floating at an air-water (left) and air-glycerol (right) inter- 
face. The diameter of the circular dish in each case is 15 cm, 
and the underlying fluid layer is 2 cm thick. 

direction || . This fracture can be observed in the partic- 
ulate scum that forms on the surface of a cup of black tea, 
which is then fractured by the addition of milk. A similar 
phenomenon is observed in pond scum when fractured by 
the ripples induced by a pebble. Here we use an interfa- 
cial particle raft as a vehicle to study the cracks induced 
by the addition of surfactant and provide a model for 
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their unusual dynamics. This study provides a window 
on the nonlinear rheology of these fluid-solid interfacial 
composites. As we will see, the dynamics of crack growth 
is governed by the rate at which surfactant can be sup- 
plied to the tip. 

Our experimental system is a densely packed mono- 
layer of non-Brownian particles || trapped at the 
interface between air and a water-glycerol mixture con- 
tained within a circular dish of diameter 15 cm and depth 
2 cm. A drop of surfactant (whose volume does not affect 
the results) is introduced somewhere in the layer with a 
needle, which when clean is benign, i.e. it does not open 
a crack. We used Polyoxyethylene Sorbitan Monooleate 
(EM Science) as the surfactant, though household deter- 
gents have the same effect. The localized reduction of 
surface tension causes a tensile stress in the particulate 
layer; a crack nucleates at the needle, grows and eventu- 
ally branches, as shown in fig. HI The composition of the 
water-glycerol mixture was varied allowing us to vary the 
viscosity of the underlying liquid, /i, from 10 -3 — 0.3 Pas. 
This variation has a pronounced effect on the speed of 
propagation of the crack (note the disparate timescales 
in fig. [1} so that the typical propagation speeds in our 
system lie in the range 10 ~ 3 — 10 _1 ms _1 . These cracks 
thus travel at speeds of between 0.2% and 20% of the 
shear wave speed, which scales as 

v 3 ~ \ —~ J 7 i~ 0.5ms" 1 , (1) 

where we have used the estimate for the shear modu- 
lus G and Poisson ratio v determined and verified in Q , 
and p was taken as the density of water and the (neu- 
trally buoyant) Pliolite. Over the large range of crack 
speeds in our experiment, we observe phenomena such 
as crack branching, kinking and the appearance of frus- 
trated side-branches, which makes our system rather un- 
usual. Indeed, they question the hypothesis that dynam- 
ical effects alone can give rise to these phenomena and 
suggest that instead heterogeneities may be the domi- 
nant cause of these effects. As an example, we consider 
the process that leads to the kinking of the crack path 
which seems to be governed by a simple phenomenolog- 
ical rule: the crack propagates in a particular direction 
through a grain of the solid until it reaches the grain 
boundary. At this point it changes direction to travel in 
whichever of the two possible cleavage directions is clos- 
est to its current direction of travel. Understanding this 
requires a consideration of the effects of disorder; indeed 
our experiments show quite clearly that crack branching 
and kinking are not due to inert ial effects alone, contrary 
to current thinking on the subject. 

The fracture of these particle rafts was observed using 
a high speed CCD camera at frame rates of up to 200s -1 . 
The images produced were subsequently analysed using 
image analysis software (DigiFlow, DL Research Part- 
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FIG. 2: Time series of the two-dimensional divergence of the 
velocity field in a raft of 90/xm particles showing the prop- 
agation of a sound wave (circular wave fronts) shortly after 
the injection of surfactant. The speed of this wave agrees 
well with that based on the estimate of the Young's modulus 
obtained by buckling experiments The orientation and 
length of the crack is represented by the solid black lines in 
each of the last three frames. 

ners) to determine the velocity field in the raft and the 
divergence of this velocity field. Typically the time series 
of the divergence field in the raft at short times (see fig. 
13 for an example) shows a wave advancing ahead of the 
crack. The typical speed of this wave is ~ 0(O.5)ms -1 , 
which is very close to the speed of shear waves in the solid 
v s , but which is considerably faster than the speed of the 
crack itself. Thus inertial effects lead to rapid dynamic 
equilibrium in the raft which learns of the elastic distur- 
bance due to the presence of the crack acoustically, even 
though the crack itself propagates much more slowly. 

Image analysis also reveals the spatial inhomogeneity 
of the response of the raft. In particular, ahead of an 
advancing crack there is a shadow zone where the veloc- 
ity of the raft is negligible, as shown in fig. 03 Since the 
speed of crack propagation is much slower than the speed 
at which the details of the crack are communicated to the 
raft and its boundaries, we assume that the displacement 
within the raft is determined instantaneously by the ge- 
ometry of the crack. Thus we may turn to the classical 
plane stress solution for the equilibrium of a plate with 
an elliptical hole subject to internal pressure 0, H^ . 
We summarize briefly here the approach followed in sec- 
tion 62 of [l3|. Since the raft is loaded in plane stress, 
the Airy stress function, (/>, satisfies V 4 (/> = 0. Letting 
z = x + iy where (x, y) is a Cartesian co-ordinate system 
centred on the hole, we introduce the complex potentials 
i/j(z) and x( z ) with </> = Re (2^(2:) + x( z ))- The bound- 
ary condition that the stress at infinity is zero leads to 
Re(i/j'(z)) = zijj"(z) + x"(z) = there. At the elliptical 



FIG. 3: The computed (left) and measured (right) velocity 
fields around the crack tip (see text). The arrows illustrate 
the instantaneous velocity at a particular point, but have been 
scaled for clarity and represent the distance that the point 
would move in the next 0.065 s. In the model (see eq. J3J), the 
geometry of the crack (solid black line on right) is represented 
by an elliptical hole of similar dimensions. The scale bar 
represents 5 mm. 
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FIG. 4: Main figure: Non-dimensional crack length as a func- 
tion of non-dimensional time for various viscosities and par- 
ticle sizes (points). Inset: Experimental confirmation of the 
scaling t~\ ~ t~l t giving the growth rate in terms of the 
different experimental parameters as in (4). 



boundary of the hole, the stress must equal the spread- 
ing pressure, S. This condition is simplified by moving 
into an elliptical co-ordinate system z = ccosh£ with 
( = £ + if] where £ is constant on a given ellipse and 
on the boundary of the hole £ = £o- This leads to the 
potentials 



S 

^0) = -(csinhC - z), xO) 



-|c 2 Ccosh2£ , (2) 



which in turn give the displacement field, (u,v), as 
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where v is the Poisson ratio (y = l/y/3 for hexagonally 
close-packed spheres interacting via a central force 
We may then obtain a theoretical estimate of the veloc- 
ity field by using the measured crack shapes at two in- 
stants of time to determine the displacement fields from 
((HJ and subtracting them. The results of this are shown 
in fig. [31 where a value of S ~ 0.3j/d has been used to 
fit theory to experiment. Integrating this fitted value of 
S over the thickness of the layer shows that the differ- 
ence in surface tensions between the crack and the solid 
is O.37. This is indistinguishable from the directly mea- 
sured difference in interfacial tensions between pure and 
contaminated interfaces. The agreement between com- 
putation and experiment shown in fig. 03 is reasonable 
given the simplicity of our model and the departure of 
the real crack from our elliptical idealisation. Over the 
points shown in fig. El the average error is 30%, while the 
angle of the shadow zone ahead of the crack tip is similar 
in both instances. 

We now move towards a quantitative description of 
the dynamics of a single crack by measuring its length, 



L c , as a function of time, t, following its initiation. It 
is observed experimentally that the cracks terminate at 
some final length, L^, which does not vary between ex- 
periments in a systematic manner. This suggests that 
the heterogeneity of packing, the initial configuration of 
the raft and the proximity of other cracks may influence 
the propagation of a single crack. We note that after 
compressing the raft to the point at which buckling sets 
in and then adding surfactant, no crack was formed: if 
there is insufficient free area within the raft for particle 
rearrangements to occur, the crack will not propagate at 
all. Conversely, once a crack has reached its equilibrium 
length, any further compression of the raft (in the di- 
rection perpendicular to the crack) leads immediately to 
the onset of buckling. In addition, the crack maintains its 
final shape for hours after the crack has ceased to prop- 
agate. This suggests that, as the crack propagates, the 
particles are consolidated to liberate liquid-gas surface 
area at the crack, while after the crack has progressed, 
the particles within the raft are in a stable jammed state 
compared to that prior to crack initiation. 

The separation of speeds between the propagation of 
the crack and the sound waves shown in fig. [21 suggests 
that the fracture itself is limited by the rate of advection 
of surfactant to the crack tip. To test this possibility, we 
used neutrally buoyant food coloring as the surfactant 
and observed that throughout the motion, the dye was 
confined to the crack area and that cracks only propa- 
gated with visible amounts of dye at the tip. This sug- 
gests that the spreading of surfactant is indeed vitally 
important for crack propagation. We now move on to 
consider the dynamics of this process. 

As the crack propagates, energy is liberated by the re- 
duction in surface tension of the crack area compared to 
the pure liquid. The dynamics of the crack propagation 
are thus governed by the balance between the rate at 
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which this energy is produced and the rate at which en- 
ergy is dissipated by the motion of the raft. There are two 
possible dissipation mechanisms. The first of these is the 
dissipation that occurs in a Blasius boundary layer, which 
is caused by the impulsive surface flow induced by differ- 
ences in surfactant concentration via Marangoni stresses. 
A time t after initiation of the crack, the typical verti- 
cal extent of this surface boundary layer is S ~ (pt/p) 1 ^ 2 
with p the dynamic viscosity of the liquid. The rate of 
viscous dissipation within the volume, V, of this bound- 
ary layer is V m ~ 2p J v (VufdV ~ p(L c /5) 2 L 2 c 5. There 
is another source of dissipation due predominantly to the 
lubrication flow in the spaces between particles: a typical 
shear rate, 7, gives a dissipation rate of ~ pL 2 dj 2 , 
an overestimate since it neglects the crucial fact that the 
liquid between particles can move out of the plane. As- 
suming that this shear rate is of the same magnitude as 
the rate of compression of the raft, we find experimentally 
that 7 ~ Is -1 so that Viub/V^u < 10 -5 <C 1, i.e. we can 
neglect the dissipation caused by inter-particle motions 
and focus on the balance between the liberation of sur- 
face energy and dissipation in the sub-surface boundary 
layer, considered by Hoult and others (jjj. Then assum- 
ing that the crack area is A c , the rate of energy liber- 
ation by the propagation of surfactant scales as AjA c , 
where A7 is the change in surface tension coefficient be- 
tween pure and contaminated interfaces. Balancing this 
with the rate of energy dissipation in the boundary layer 
V hU ~ p(L c /^/pt/p) 2 A c ^/pt/p, yields 



A7 2 
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,3/4 



(4) 



In fig. 21 we show that the scaled crack length L = 
Lc/Loq ~ t 3 / 4 over more than a decade in the scaled 
time r = t/ti/2 where ti/ 2 is the time taken for the crack 
to grow to half of its final length. Eventually, the crack 
slows down as it propagates into regions of densely packed 
particles and finally stops. In the inset to fig. 0] we also 
see that t~j 2 ~ t~^ t , where t nat = (/ipL^ c /A7 2 ) 1//3 is the 
natural time-scale from Q. We observe that the bound- 
ary layer extends through the depth, h, of the fluid sub- 
layer when t > h 2 p/p. For the raft shown in the right- 
hand panel of fig. p = 1.4 Pas so that this time is 
0(1 s), explaining why this crack propagates even slower 
than we would expect on the basis of Q. 

Our experiments on the surfactant-induced failure in 
an interfacial particulate raft show some unusual features 
such as crack branching and kinking even though the 
fractures are limited by the drag induced by the forma- 
tion of a viscous boundary layer in the underlying fluid, 
rather than solid inertia. They also confirm the solid- 
like response of densely packed particles trapped at a 
liquid-gas interface. Understanding the means by which 
fractures occur and the dynamics of their propagation in 
these systems may have very real benefits in, for example, 



stabilising foams, drops and bubbles by the addition of 
particles. One particularly interesting possibility lies in 
the delivery of drugs by inhalation, where trojan parcels 
of drugs in drops are coated with particulate rafts 
Our work suggests that the delivery process might be 
expedited by using the surfactant lining of the lungs to 
induce the rupture of the particle shell. On a different 
note, this transport-induced weakening and growth of a 
crack is similar to what occurs in corrosion cracking and 
hydrogen embrittlement in solids, and thus our system 
may also provide a two-dimensional model to study some 
of these complex phenomena. 
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